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Abstract 

In the standard seesaw mechanism, finite corrections to the neu- 
trino mass matrix arise from 1-loop self-energy diagrams mediated by a 
heavy neutrino. We study in detail these corrections and demonstrate 
that they can be very significant, exceeding in several cases the tree- 
level result. We consider the normal and inverted hierarchy spectra 
for light neutrinos and compute the finite corrections to the different 
elements of the neutrino mass matrix. Special attention is paid to 
their dependence with the parameters of the seesaw model. Among 
the cases in which the corrections can be large, we identify the fine- 
tuned models considered previously in the literature, where a strong 
cancellation between the different parameters is required to achieve 
compatibility with the experimental data. As a particular example, 
we also analyze how these corrections modify the tribimaximal mixing 
pattern and find that the deviations may be sizable, in particular for 
#13. Finally, we emphasize that due to their large size, the finite cor- 
rections to neutrino masses have to be taken into account if one wants 
to properly scan the parameter space of seesaw models. 
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1 Introduction 



Neutrino oscillation experiments have firmly established that neutrino have 
tiny but non-zero masses and that the mixing in the leptonic sector is in 
sharp contrast with the small mixing that characterizes the quark sector 
[TJE]. From a theoretical perspective the smallness of neutrino masses can be 
well understood within the seesaw model [3H7], in which the fermion sector 
of the Standard Model is extended by adding new electroweak fermionic 
singlets (standard seesaw model). In this framework, light neutrino masses 
are generated via mixing with the singlet states and their smallness can 
naturally be explained if the singlets masses are very large. 

The determination of the regions of parameter space consistent with low 
energy neutrino observables in the seesaw model typically relies on parametriza- 
tions of the neutrino Yukawa couplings [8,9]. Once the Yukawas are properly 
parametrized, such regions are found by doing numerical scans in which the 
neutrino experimental data is used as an input. This procedure is always 
based on the tree-level light neutrino mass matrix and fails if in some regions 
of the parameter space the one-loop corrections to the tree-level mass matrix 
turn out to have sizable values. 

The one-loop corrections to the seesaw light neutrino mass matrix were 
first discussed in |1U| in a general setup with an arbitrary number of singlets, 
lepton doublets and Higgs doublets. They were later analyzed in ref. [UJ in 
a particular realization in which, due to a particular Yukawa mass matrix, 
light neutrino masses vanish at tree-level and are entirely generated by the 
one-loop corrections. Subsequently, the renormalization of general theories 
with Dirac and/or Majorana neutrinos was carried out in |12] and additional 
studies were done in references [T5HT7] . 

Loop corrections are of two types: renormalizable and intrinsically finite. 
The renormalizable pieces consist of corrections to the tree level parameters 
already present in the seesaw Lagrangian, and are suppressed with respect 
to the tree level piece by extra Yukawa couplings and by the loop factor 
1/167T 2 . The finite parts instead are corrections to the vanishing elements of 
the tree-level mass matrix for the neutral fermions, and are only suppressed 
by the loop factor. Thus, they are potentially large. 

The aim of this paper is to quantify the impact that the finite one-loop 
corrections might have on the effective light neutrino mass matrix. We con- 
sider the most general standard seesaw model and numerically analyze the 
importance of these corrections for the different mass matrix elements as well 
as for the neutrino mass eigenvalues and mixing angles, differentiating in our 
discussion between the normal and the inverted light neutrino mass spectrum. 
We will show that the corrections can range over several orders of magnitude, 
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depending on whether one relies or not on models where consistency with the 
measured neutrino masses and mixing angles requires strong cancellations in 
the tree-level neutrino mass matrix. Indeed, as we will discuss, once one-loop 
corrections are taken into account these models are barely reconcilable with 
data. Barring these cases, we will prove that the finite corrections are usually 
of order 20 — 40% but they may exceed the tree- level value for certain entries 
that can be strongly suppressed. In order to make reliable predictions in the 
seesaw, therefore, the finite one-loop corrections should be included. 

The rest of the paper is organized as follows: in section [2] we define our 
notation and briefly describe the seesaw model at tree-level, including its 
standard parametrization. We discuss the finite 1-loop corrections in sec- 
tion [21 In sections H] and our main results are presented -the calculation of 
the finite corrections for the normal and inverted light neutrino mass spec- 
trum. Then, we consider the particular case of tribimaximal mixing in section 
[6]and we determined how such mixing pattern is modified by the 1-loop finite 
corrections. A brief discussion of our results and some comments on their 
possible phenomenological implications is given in section [71 In section [HI 
we draw our conclusions and summarize our findings. For completeness, in 
the appendix we present the details of the calculation of the finite one-loop 
corrections. 

2 The standard seesaw model at tree level 

In the standard seesaw model, three fermionic electroweak singlets Nj^ are 
added to the Standard Model. In the basis in which the matrix of charged 
lepton Yukawa couplings and the singlet mass matrix are diagonal the La- 
grangian accounting for the new interactions can be written as 



where (f) T = is the Higgs electroweak doublet, £l are the leptonic 

SU(2) doublets, C is the charge conjugation operator and A is a Yukawa 
matrix in flavor space. In the left-handed chiral basis iil T = (i/l, (Nr) c ) t , 
and once electroweak symmetry breaking is taken into account, the neutral 
fermion mass terms can be written as 





- C F o = -tl^CM n L + h.c. 

where Ai, the 6x6 neutral fermion mass matrix, is given by 



(2) 




(3) 
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with M D = vX (v = V2M w /g ~ 174 GeV). The mass spectrum is obtained 
by rotating the fields to the mass eigenstate basis, denoted by x%i y i a the 
unitary matrix U: 

XL = UW . (4) 
In this basis, the Lagrangian mass terms, equation (J2J), become 

-jC F o = ~xMP LX + h.c. (5) 

where 

U T M U = M = diag(m Xl , . . . m X6 ) (6) 

and X i are the physical Majorana neutrino fields. Note that by decomposing 
the matrix U as ITUlfTHlfTT] 

U L 

u 



U=( T ^) (7) 



the and states can be expressed as 

v Li = U Lij P LXj , N Ri = U Rij P l0 C j . (8) 

In the seesaw limit (Md <C Mr), the diagonalization of the mass matrix ([3]) 
gives rise to a split spectrum consisting of three heavy states with masses 
Mr. and three light states with an effective mass matrix, m^ ree) , given by 

m (tree) = _ M T Md _ (g) 

The light neutrino mass spectrum, mixing angles and CP violating phases 
-the so-called low-energy observables- are obtained from this matrix after 
diagonalization: 

Ujm^U^m,, (10) 
where is the leptonic mixing matrix parametrized according to 

U, = U^ 23 )U^i3,5)U,(M x diag(e-^,e^ 2 ,l) (11) 

with 5, (px 2 being respectively the Dirac and Majorana CP violating phases 
and Ue(6) rotation matrices. 

The determination of the seesaw parameters compatible with neutrino 
experimental data relies on parametrizations of the Yukawa couplings or, 
equivalently, of the Dirac neutrino mass matrix, Md = v X. In the numerical 
analysis of the finite one-loop corrections, we have used the most common 
parametrization of the seesaw, the Casas-Ibarra parametrization [8J. In this 
parametrization, the most general Md compatible with eq. (jH]) is given by 

M D = iM^Rmy 2 U/ , (12) 
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where R is any orthogonal matrix. This matrix can be written as a rotation 
matrix determined by three complex angles. 

The neutrino mass eigenvalues and the mixing matrix entering into this 
equation are strongly constrained by experimental data whereas the masses 
of the singlet neutrinos and the matrix R are entirely free parameters of the 
seesaw model. In the numerical treatment of our results (sections H] and \5§ 
we also impose the perturbativity condition suggested recently in |18| : 

Tr [A + A] < 3. (13) 

Now that we have reviewed the seesaw mechanism at tree level, let us take a 
look at the 1-loop corrections to neutrino masses. 

3 Finite one loop corrections to the neutral 
fermion mass matrix 

In the standard seesaw, the one-loop corrections to the v — N mass matrix 
are determined by the neutrino interactions with the Z boson, the neutral 
Goldstone bosons (G°), and the Higgs boson (h°) -see appendix |A1 All 
together, in addition to the correction involving the standard model leptonic 
charged current, they define the one-loop two point function — iS(p) |17| . 

Once the one-loop corrections are taken into account the neutral fermion 
mass matrix is given by 

M = M (tIcc) + M (1 - ,oop \ (14) 
where the 1-loop contribution can be decomposed as 

Notice that the O3X3 matrix appearing at tree-level is replaced by the contri- 
bution <5Ml, which among all the sub-matrices in J\A (1Aoop) is the dominant 
one [17J. 

Neglecting the subdominant pieces in 7W (1 " loop) and after block diagonal- 
ization of the neutral fermion mass matrix, the effective light neutrino mass 
matrix, up to one-loop order, can be written as 

= m|T e) + m^ loop) = -M£ Mr 1 M d + 5M L . (16) 

The sub-matrix 5Ml and all the other sub-matrices entering in 7Vl (1 " loop) are 
entirely determined by the self-energy functions S£(p 2 ) (see appendix |AJ via 
the diagonalization relation ()6]): 

^d-ioop) = U *I]S(p 2 )U t . (17) 
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Accordingly, the finite contribution is given by 



<5M L = U L *S^(/)U^ = U L *E£(0)U L t , (18) 

where we have used the fact that S£ can be evaluated at zero external 
momentum [IT]. The self-energy functions S£(0) are determined by three 
Feynman self-energy diagrams involving the Z, the neutral Goldstone boson 
G° and the Higgs boson h°. Each diagram contains a divergent piece but 
when summing up the three contributions the result turns out to be finite, 
as it has to be since there are no counterterms that would allow to absorb a 
possible divergence (see appendix |2] for more details). The final expression 
for the finite one-loop correction is given by [T7] 



5M T , = MlM R x I -A 



64tt 2 M^ 



M 



D 



(19) 

Notice that this correction is not suppressed, with respect to tree- level result, 
by additional factors of Md/Mr. Thus, it is expected to be smaller than the 
tree-level mass term solely by a factor of order (167r 2 ) -1 ]h(Mr/Mz)- 

In spite of the similar structure of the 1-loop correction and the tree-level 
result, they are not proportional to each other unless the heavy neutrinos are 
degenerate Mr oc I. Hence, one could in principle have that m[ / tloc) = and 
that neutrino masses arise entirely from 1-loop effects, as proposed in [llj. 
Such models, however, are rather contrived and will not be discussed in the 
following. We are interested, instead, in the generic modifications to the 
neutrino mass matrix induced by the 1-loop corrections. 

To evaluate these corrections, we first find sets of and Mr compatible 
with the experimental data at tree-level -using equation (TT2]) - and then use 
them to evaluate Specifically, we generate the diagonal matrix of light 

neutrino masses (according to the desired spectra: normal or inverted) and 
the mixing matrix such that they are compatible with neutrino data. For 
simplicity, the phases in were assumed to vanish. Then, we randomly 
generate the three masses of the heavy states (in the range 1 TeV to 10 12 
GeV) and the elements of the orthogonal matrix R. From equation ()12|) . we 
can then obtain Md, which together with the generated Mr allows us to 
evaluate 

smS- The size of the corrections is then determined by the ratio 
between the contributions up to 1-loop order and the tree-level result for the 
different elements of the neutrino mass matrix, (m^ tree) + m[ / 1 " loop) ) /m^ tree) . 



Alternatively, one could choose the seesaw parameters so that at 1-loop they are 
compatible with the experimental data. Both procedures give rise to the same effects. 
b In our analysis we fix myp = 150 GeV. 
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In the next two sections our main results are presented: we compute the 
corrections to the neutrino mass matrix in the seesaw model for the two 
different kinds of light neutrino spectra, with normal and inverted hierarchy. 



Corrections for the Normal Hierarchy spec- 
trum 



If the spectrum of light neutrinos has a normal hierarchy (m V3 = a/ Am^ tm , 
m V2 = a/ Am 2 oi , m Ul <C m U2 ,m u ^), the elements of the neutrino mass matrix 
take values within the following ranges: 

/ (2.5,5.5) x 1(T 12 (2.3,9.8) x 1(T 12 (-3.3,4.9) x KT 12 \ 
rrC P = - (2.0,3.4) x 10" 11 (1.9, 2.3) x 10 -11 GeV 

\ - - (2.1,3.4) x 10- 11 / 

(20) 

as the oscillation parameters vary within their 2-cr experimentally allowed 
intervals [HE]. Since the matrix is symmetric, we only show the six inde- 
pendent matrix elements. Notice, in particular, that the element (1,3) is 
the only one that can vanish in this case. It can be easily checked that this 
can happen if # 13 is between 4° and 6°. All other elements vary within a 
relatively small range -not so small for (1,2)- between 10~ n and 10~ 12 GeV. 
Since the corrections to the neutrino mass matrix are not proportional to the 
matrix element itself, the correction to the element (1, 3) could easily exceed 
its tree- level value. 

With the aim of facilitating the study of these corrections and the under- 
standing of their origin, we will divide our analysis in two parts depending 
on what is assumed for the orthogonal matrix R. First it is taken to be real 
and then the most general case is considered, a complex matrix. The number 
of free parameters will therefore increase as we move from the first case to 
the second. 



4.1 R real 

If R is real, the parameters needed to evaluate the 1-loop correction to the 
neutrino mass matrix are the oscillation parameters, the three masses of the 
right handed neutrinos (Mi), and the three angles that parametrize R. To 
obtain the numerical results below, we vary the neutrino mixing angles and 
mass squared differences within their 2a ranges, and we randomly choose Mi 
between 1 TeV and 10 12 GeV and the angles in R between and 2n. 
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Figure 1: The ratio between the 1-loop and the tree-level result for the different 
elements of the neutrino mass matrix as a function of Mi . It has been assumed 
that R is real and that light neutrinos have a NH spectrum. 



The resulting corrections to the matrix elements are shown in figure [T] as 
a function of the mass of the lightest heavy neutrino, M±. We see that they 
are similar for the elements (1,1), (2,2), (2,3), and (3,3), increasing with 
Mi and reaching values up to order 30%. Those for the element (1,2) are 
slightly different, reaching values as large as 40% or 50% as well as —20%. 

The corrections to the element (1,3), on the other hand, can be quite 
large for a significant fraction of models. As anticipated, this result is due 
to the fact that the element (1,3) can be very small so it may receive a 
huge fractional correction. It must be noticed in that case, however, that a 
large value of (m£ ree) + m^~ loop) )/m^ ree) does not necessarily imply a signifi- 
cant deviation in the expected value of the neutrino observables -the mass 
eigenvalues and the mixing angles. For that reason it is important to study 
the effect of the corrections on both the matrix elements and the predicted 
observables. We will do so in the next section, where we consider the most 
general case: R complex. 

Notice then that even in the case R real, where no large parameters are 
introduced, the corrections to neutrino masses can be quite important. If 
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M, [GeV] Mj [GeV] Mj [GeV] 



Figure 2: The ratio between the 1-loop and the tree-level result for the different 
elements of the neutrino mass matrix as a function of Mi . It has been assumed 
that R is complex and that light neutrinos have a NH spectrum. 



Mi > 10 9 GeV they are expected to be larger than about 15% and they 
could easily reach 25% or 30%. 

4.2 R complex 

This case is not only the most general one, but it is also well motivated by 
leptogenesis. In fact, in the case of unflavored leptogenesis, the phases in R 
are the ones responsible for the CP-asymmetry in the decays of the heavy 
neutrinos and ultimately for the generation of the baryon asymmetry. 

In this case, the three angles parametrizing R are complex numbers, 
with a certain magnitude and a given phase. In our analysis, we allow these 
complex angles to have an arbitrary phase and we restrict their magnitude 
to be smaller than 30. Since cosh 3 ~ sinh3 ~ 10, R can have elements 
at most of order 10 3 . When the elements of R are significantly larger than 
1, \Rij\ ^> 1, one obtains the so-called fine-tuned models. In them, strong 
cancellations between the different terms in equation are required to 

c One can relax this assumption to obtain even larger effects. 
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10 100 1 10 100 '1 10 

IR.J IR..I IR..I 



lj max y max lj max 

Figure 3: The ratio between the 1-loop and the tree-level result for the different 
elements of the neutrino mass matrix as a function of the largest element of the R 
matrix. It has been assumed that R is complex and that light neutrinos have a NH 
spectrum. 



obtain compatibility with the experimental data. Since the corrections to 
the neutrino mass matrix depend on R, the loop suppression in (|T9|) can 
be easily overcome by the large elements in R, yielding a correction that is 
significantly larger than the tree level result. 

Figure [2] shows the corrections to the different matrix elements for R 
complex. We see that, in fact, the 1-loop contribution can exceed, for all the 
matrix elements, the three-level result by several orders of magnitude. In 
that case, there is no doubt that the corrections will have a huge impact on 
the predicted neutrino mass eigenvalues and mixing angles. 

It is indeed the large elements present in R that make possible a 1-loop 
correction much larger than the tree- level result. We illustrate this fact in 
figure [31 which shows the size of the corrections as a function of the largest 
element of the R matrix, \Ra\ . Notice that when this element is of order 

" I J I max 

1 the corrections are usually small (those to the element (1,3) being the 
exception) but they increase with it reaching two orders of magnitude or 
more for IFLA around 100. With such huge corrections, the agreement 

I L J I max ° ' ° 
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Figure 4: The neutrino mixing angles at 1-loop as a function of the largest element 
of the R matrix. It has been assumed that R is complex and that light neutrinos 
have a NH spectrum. The region between the two dashed (red) lines is consistent 
with current experimental data at 2a. 



between the tree-level seesaw formula and the neutrino data assumed in the 
parametrization becomes meaningless. In fact, as illustrated in figures H] and 
the oscillation parameters may deviate significantly from their observed 
values once the 1-loop corrections are taken into account. 

Figure H] displays the 1-loop mixing angles, those obtained from the diag- 
onalization of the neutrino mass matrix at 1-loop, as a function of \Rn\ 

" ' I J I max 

The region consistent with the experimental data at 2cr is the area between 
the two dashed lines. Notice that all the angles, which were chosen to be 
consistent with the data at tree level, can at 1-loop become much larger than 
allowed by present observations. 

Similarly, we see in figure [5] that the neutrino mass squared differences 
at 1-loop can vary over several orders of magnitude. A fact that is in clear 
contradiction with current experiments. 

As we have seen, for the neutrino spectrum with normal hierarchy the 
corrections to the neutrino mass matrix can be quite important. The ma- 
trix element (1, 3), in particular, can receive very large fractional corrections 
independently of R if m ul3 is suppressed. In addition, all matrix elements 
as well as the neutrino mass eigenvalues and mixing angles are expected to 
receive significant corrections when the elements of the R matrix are larger 
than one. In such case, the inclusion of the 1-loop corrections is mandatory. 

Next, we analyze the importance of these corrections for the neutrino 
spectrum with inverted hierarchy. 
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Figure 5: The neutrino mass squared differences at 1-loop as a function of the 
largest element of the R matrix. It has been assumed that R is complex and that 
light neutrinos have a NH spectrum. 



5 Corrections for the Inverted Hierarchy spec- 
trum 

If the spectrum of light neutrinos has an inverted hierarchy (IH), m U2 = 
a/| Am 2 im |, m ui = ^/\/\rn? atm \ — Am^, m U3 <C m V2 , m Ul , the elements of the 
neutrino mass matrix take values within the following ranges: 

(4.5,5.1) x 10- 11 (-8.5,1.7) x 10- 12 (-8.5, 1.3) x 10" 12 \ 

(1.7,3.3) x 10- 11 (-2.5,-2.0) x 10- 11 GeV 

(1.8,3.4) x 1CT 11 / 
(21) 

as the oscillation parameters vary within their 2-a experimentally allowed 
intervals [U|2]. Notice that in this case the elements (1, 2) and (1, 3) are both 
allowed to vanish, an event that can happen if 6*13 is smaller than about 2°. 
We expect, therefore, large fractional corrections to the entries (1,2) and 
(1,3) of the neutrino mass matrix independently of R. As before, we will 
divide the analysis of the 1-loop corrections into two parts: R real and R 
complex. 

5.1 R real 

The finite corrections for R real are shown, as a function of M 1; in figure 
O The range of variation is approximately the same for the elements (2, 2), 
(2, 3), (3, 3), and (1, 1), reaching maximum values of order 30% independently 
of M\. The minimum value of the correction, on the other hand, is seen to 
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Figure 6: The ratio between the 1-loop and the tree-level result for the different 
elements of the neutrino mass matrix as a function of Mi . It has been assumed 
that R is real and that light neutrinos have a IH spectrum. 



increase with Mi. The entries (1,2) and (1,3) may feature large fractional 
corrections, a consequence of the vanishing matrix elements at tree-level. 
Comparing these results with those obtained in the previous section, it is 
evident that the type of light neutrino spectrum does not have a decisive 
impact on the generic size or behavior of the finite corrections. 

5.2 R complex 

In the most general case of a complex R matrix, the corrections tend to be 
quite large for all entries, as illustrated by figure They can easily reach 2 
or 3 orders of magnitude above the tree-level value, being typically larger for 
the (1,2) and (1,3) matrix elements. They may also give rise to cancellations 
between the tree-level and the 1-loop contribution, such that the full result 
at 1-loop could only be a small fraction of the tree-level result -see e.g. the 
points around 0.1 in the figure. 

These large fractional corrections to the elements of the neutrino mass 
matrix translate into important deviations in the neutrino mass eigenval- 
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M, [GeV] M, [GeV] M t [GeV] 



Figure 7: The ratio between the 1-loop and the tree-level result for the different 
elements of the neutrino mass matrix as a function of Mi . It has been assumed 
that R is complex and that light neutrinos have a IH spectrum. 



ues and the neutrino mixing angles, just as for the spectrum with normal 
hierarchy -see previous section. 

For a light neutrino spectrum with inverted hierarchy, therefore, the cor- 
rections are even more important than for the normal hierarchy spectrum, 
as there are two different matrix elements that can receive large fractional 
corrections independently of R. If R contains large numbers then the correc- 
tions to all matrix elements are usually significant for both the normal and 
the inverted hierarchy spectrum. 

6 A specific example: tribimaximal mixing 

We would like now to apply the ideas discussed in the previous sections to 
a particular and well-motivated scenario: seesaw models with tribimaximal 
mixing (see e.g. [19|). In scenarios with tribimaximal mixing, the neutrino 
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10000 le+06 le+08 le+10 le+12 

Mj [GeV] 



Figure 8: The mixing angle 623 at 1-loop as a function of the lightest heavy neu- 
trino. At tree level #23 was assumed to be 45°, in agreement with the tribimaximal 
mixing pattern. 



mixing matrix is given at tree level by 

/ v^Ts 1/V3 \ 

U t = -1/VE 1/VS -1/V2 . (22) 
V -1/V6 l/y/3 l/y/2 ) 

As a result, the mixing angles take the values 6 12 = 35.3°, #23 = 45°, 9% 3 = 
at tree level. One may therefore wonder how they would change once the 
finite one-loop corrections to the seesaw neutrino mass matrix that we have 
studied are taken into account For simplicity, we will limit ourselves in this 
section to the normal hierarchy spectrum and to the case R real. Larger 
corrections are expected if R is complex. 

Figure [S] shows the 1-loop value of 623 as a function of the lightest heavy 
neutrino mass. We see that it can deviate from its tree-level value by up to 
2.5 degrees and that the maximum deviation decreases with M%. In fact, for 

d Additional corrections from other sources may also be important but are not consid- 
ered here. 



15 




10000 le+06 le+08 le+10 le+12 

Mj [GeV] 

Figure 9: The mixing angle 612 at 1-loop as a function of the lightest heavy neu- 
trino. At tree level Qyi was assumed to be 35.3°, in agreement with the tribimaximal 
mixing pattern. 



M\ > 10 12 GeV the correction is smaller than half a degree. This dependence 
with Mi is very different to that observed for the matrix elements but it is 
consistent with it. Since the mixing angles are determined by ratios between 
different matrix elements, in order to have a sizable variation in the mixing 
angles, the structure of the neutrino mass matrix should vary significantly at 
1-loop, implying that the correction should no be proportional to the identity. 
Hence, in models where there is a large hierarchy between the masses of the 
heavy neutrinos, the corrections to the mixing angles are larger. That is 
exactly what is observed in figure [BJ 

The 1-loop corrected value of 9 12 is shown in figure [9] as a function of M\. 
The variation in this case is smaller and it also decreases with M\. More 
interesting for the phenomenology of neutrinos and for future experiments is 
the correction to # 13 , which is exactly zero at tree level. Figure [TU] shows that 
the 1-loop corrected value of 9\$ can reach almost 2 degrees, corresponding 
to sin 2 #13 ~ 10~ 3 , for Mi around 1 TeV. As expected, the maximum correc- 
tion decreases with Mi, amounting to about 1 degree (sin 2 #13 ~ 3 x 10~ 4 ) 
for Mi ~ 10 8 GeV. Given that a neutrino factory could be sensitive to 
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10000 le+06 le+08 le+10 le+12 

Mj [GeV] 



Figure 10: The mixing angle #13 at 1-loop as a function of the lightest heavy neu- 
trino. At tree level #13 was assumed to vanish, in agreement with the tribimaximal 
mixing pattern. 

sin 2 6*i3 ~ 10~ 5 [20J, these corrections are certainly within the reach of future 
experiments. 

7 Discussion 

The huge corrections we have found for fine-tuned models are not surpris- 
ing. It is well-known in the literature that radiative corrections to fine-tuned 
models should spoil the tuning imposed at tree-level between the different 
parameters (see e.g. |21] ) and that it is unnatural to expect otherwise. We 
have explicitly shown that that is the case. The 1-loop corrections to the 
neutrino mass matrix in fine-tuned seesaw models are so large that the com- 
patibility between the tree-level seesaw formula and the experimental data 
becomes irrelevant. The inclusion of the 1-loop corrections is in such case 
necessary. 

As a way of avoiding such large corrections, one may think of including the 
1-loop correction into the seesaw parametrization from the very beginning. 
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That is, one would like to look for the most general solution to equation 
([To 7 ]) rather than to equation Doing so, however, would not completely 
solve the problem. We would, instead, be constructing fine-tuned models at 
the 1-loop level, which are expected to receive large corrections at 2-loops. 
Fine-tuned models, it seems, are better avoided. 

One possibility to do so is simply to restrict from the very beginning 
the magnitude of the complex angles that parametrize R. If they are such 
that < 1 then no fine-tuning occurs and the corrections are usually 

under control. This additional restriction, however, has not been taken into 
account in previous analysis. And it was recently suggested in [18J that 
a fair scan of the seesaw parameter space is one in which no restriction 
beyond the perturbativity of the neutrino Yukawa couplings, which we have 
implemented, is imposed. Our results clearly demonstrate, on the contrary, 
that perturbativity is not enough to guarantee the stability of the neutrino 
mass matrix under radiative corrections and that wrong results can easily 
be reached if the 1-loop corrections are not taken into account. In fact, as 
we have seen, a significant fraction of models which are compatible with the 
data at tree level are no longer so once the 1-loop corrections are considered. 
Thus, when one is randomly scanning the parameter space of the seesaw 
model, it is necessary to include the finite corrections to neutrino masses. 

Another possible way out is the use of a different parametrization. An 
alternative to the R parametrization that has been used in previous works 
is the Vl parametrization. In it, is written as 

M D = V R M D V L , (23) 

where Md is a diagonal matrix defined by the eigenvalues of Md (real and 
positive) and Vl,r are unitary matrices determined by three rotation angles 
and three complex phases. In this parametrization, Md, Vl and the neutrino 
data are used as inputs. Using equation f )23|) and the effective light neutrino 
mass matrix the following relation is obtained 

M^V'mrvtM, 1 = V^ 1 V R , (24) 

which allows us to determine Vr and Mr for a given set of input parameters. 
We have computed the corrections to the neutrino mass matrix also in this 
parametrization, and have observed that the results for R = I and R real are 
easily reproduced for Vl = I and Vl real. In particular, the large corrections 
for certain matrix elements are obtained there too. An important difference 
occurs, however, for R complex. Due to the different way in which the Vl 
parametrization samples the parameter space of the seesaw model, it is way 
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more difficult to find fine-tuned models, with the consequence that models 
with very large corrections are rather scarce in the Vl parametrization. 

In a future publication, we will discuss additional implications of these 
corrections, including their evaluation in supersymmetric scenarios as well as 
their possible effects in leptogenesis and lepton flavor violating processes. 

8 Conclusions 

The seesaw model is one of the most appealing extensions of the Standard 
Model that can explain neutrino masses. In this model, the mass matrix of 
light neutrinos receives finite corrections from 1-loop diagrams mediated by 
the heavy neutrinos. We considered the two different kinds of light neutrino 
spectra, hierarchical and inverted, and computed the corrections to the en- 
tries of the neutrino mass matrix as a function of the seesaw parameters. We 
found these corrections to be quite important, exceeding in several cases of 
interest the tree level result by orders of magnitude. Two different reasons 
were identified as leading to a large correction: an unusually suppressed tree- 
level result and an R matrix with elements much larger than 1. Examples 
of the first case are the corrections to the matrix element (1,3) for NH neu- 
trinos and to (1,2) and (1,3) for IH neutrinos. The second case can occur 
for R complex and includes the so-called fine-tuned models considered in the 
literature. Since these corrections can be large, models that at tree-level are 
compatible with the experimental neutrino data will not necessarily be so 
at the 1-loop level, modifying in a significant way the viable regions in the 
parameter space of the seesaw model. As a particular example, we studied 
the corrections to the mixing angles in seesaw scenarios with tribimaximal 
mixing and show them to lead to observable effects in future experiments. 
We stressed, therefore, that because of their size and importance, these cor- 
rections must necessarily be taken into account in the study and analysis of 
seesaw models. 
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Figure 11: Self-energy diagrams accounting for <5Ml 

A Finite self-energy functions 

In this appendix we present the calculation of the finite 1-loop corrections 
<5Ml discussed in section [3J we will closely follow ref. [T7]. The self-energy 
function El(0), that determines these corrections, can be written as 

- z£S(0) = -i [s* (z) (0) + SE (GO) (0) + £f 0) (0)] , (25) 

where the ' ' (0) functions arise from the self-energy Feynman dia- 
grams (evaluated at zero external momentum) involving the Z, the neutral 
Goldstone boson G° and the Higgs h° shown in fig. [TTJ The calculation of 
these functions is determined by the coupling of the Z with the Majorana 
eigenstates \: 

C z = y~Z,xY [PlCUl+Ul) - Pk(U l t U l *)] x , (26) 

(c w = cos8 w with 9 W the weak mixing angle), the couplings with the Higgs 
boson, derived from the Lagrangian (JTJ, 

- C h o = ^=h° X [OiP L + O s R P R ] x , (27) 

where the couplings OfJ r are given by 

0!=i4AU l + U£A t U^ (28) 
Ol =U^AtU R + Ul\*Vi , (29) 

and finally the couplings with G° that can be obtained from the Lagrangian 
(T27|) by replacing r — > —iO^ R . From these Lagrangians and the diago- 

nalization relation f)17p it can be seen that S| (0) contributes only to <5Ml 

q(G° h°) 

whereas Sl ' (0) contribute to all the block matrices of the 6x6 neutral 
fermion mass matrix. The contributions of these self-energies to <5Ml can be 
identified by means of the relation (I18p . 
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Using dimensional regularization (d — 4 — e) and working in the gauge 
the Z self-energy function is found to be 



S£ (Z) (0) = U L T [5M L ^(1) + <5M L ^(2,1) + 5M L ^(2,2)] U L , (30) 

where the different matrices can be expressed in terms of the Passarino- 
Veltman function i? (0, m 2 ,, m 2 ) [22], namely 

SM^ z \l) = -^-f(4 - e)U L * M B (0, Ml, M 2 ) U L + , (31) 



'w 

-.2 



<fM L < z >(2, 1) = - 2 % A/f2 U L * M 3 B (0,Z z M 2 z ,M 2 ) iV , (32) 

<fM L < z >(2,2) = - f U L * .M 3 B (0, Mf, .M 2 ) . (33) 
As regards the G° and /i° self-energies they are given by 

£SW = Ul^Ml^Ul (X = G°, h°) (34) 
with 5M L ( - X ' ) given by 

5M l (G0) = ^ O ^2 U l* A4 3 B (0,^M1, A4 2 ) U L t (35) 

In the calculation of the above expressions we have used the relation 

U R f M d = M U L f , (37) 

that follows from the diagonalization relation fl6]) and the unitarity con- 
straints of the matrix U. 

Some words are in order regarding these results. Corrections ( |32|) and ( 135]) 
cancel, ensuring the gauge invariance of the result. The Passarino-Veltman 
function B has a finite and infinite part]f|, the infinite piece in 5Ml^(1) 
cancels due to the constraint 

U L * M U L = , (38) 

c Here by infinite part we mean B^ 11 ^ (0, nip, mf ) = 2e~ 1 — 7 + 4tt + 1. 
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whereas the divergent pieces in <5M L (Z) (2,2) and 5M L {h0) cancel amone; 
thenf], thus demonstrating that #Ml is finite as anticipated in sec. [21 Taking 
into account that the finite part of Bq can be recasted as 



B{{0,m 2 ,ml) 



1 1 ( m W , 2 

2, 2 7 m — + lnm l 

mymQ — 1 \TOq / 
rn i/ m o m 



+ lnm 2 . 



the finite parts of 5M L {z \l) and 5M L (Z - , (2, 2) combine to yield 



(z). 



(39) 



64tt 2 M2 I M? 



-II lnl^]U L t. (40) 



Finally the finite contribution from the Higgs self-energy function reads 



5M L ^ = - t— u L *M 



-i 



64tt 2 M 



111 



U T t 



(41) 



The finite correction 5Ml, discussed in sec. [21 is obtained from the 
dominant parts of eqs. ( 14"0"]) and (|4ip (order M^ 1 ). These pieces can be 
extracted by using eq. (1371) an d by taking into account that in the seesaw 
limit Md <C Mr, in the basis for which Mr is diagonal, the matrix Ur can 
be written as 



U 



R 



-e f u,, i) 



(42) 



where f = M D T M R 1 . 
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